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Abstract 

We calculate the shear viscosity r/ of a quark plasma through the two-flavor Nambu-Jona- 
Lasinio (NJL) model in the chiral limit, at finite temperature T, and baryon number chemical 
potential /t_b. We solve the Boltzmann equation by using in the collision term cross sections that 
are correct to leading order in the coupling constant. We find the ratio of shear viscosity to entropy 
density above and slightly below the chiral transition temperature to vary from 1.5 to 13 times 
the conjectured lower bound of (47r) _1 depending on the chemical potential to temperature ratio. 
Ratio of shear viscosity to entropy density is found to be a monotonically increasing function of 
the ratio of chemical potential to temperature for temperatures between 185 and 300 MeV. 


1 Introduction 

Asymptotic freedom property of Quantum Chromodynamics (QCD) allows rigorous perturbative cal¬ 
culations of transport coefficients to be done at temperatures high compared to the QCD scale A qcd- 
Recent work over the last two decades has resulted in the calculation of transport coefficients to leading 
order in the strong coupling constant for high temperature QCD [U El la g. The drawback of these 
perturbative methods is that they require the strong coupling constant to be perturbatively small. 
However, at the temperature of 200 MeV attained at the Relativistic Heavy Ion Collider (RHIC) or 
the higher temperature of 490 MeV that is expected to be reached at the Linear Hadron Collider 
(LHC) this condition is not fulfilled. This motivates the search for alternative methods to calculate 
these transport coefficients at the relevant temperatures. 

One such alternative method is the use of models that possess symmetries of low energy QCD to 
calculate the transport coefficients of interest. For two flavors of massless quarks these symmetries 
are SU(2) V x SU(2) A x U(1) B . The proposed model, of course, should also be able to account for 
the spontaneous breaking of the axial symmetries, and the invariance of the ground state under the 
diagonal subgroup SU(2) V x U(1) B . Corresponding to each broken (axial) generator a Goldstone boson 
appears, which are then naturally identified with the neutral and charged pions. 

NJL model [5] which was proposed quite a time ago is one such model. It gives a realistic description 
of low energy QCD. It exhibits a chiral phase transition through which the quarks (and anti-quarks) 
acquire masses as the temperature is lowered below the critical temperatur(0. Even though it was 
put forth as a theory of nucleons, it was reinterpreted as a theory of quarks, and extended to finite 
temperature and nonvanishing baryon number chemical potential. 

1 Critical temperature will always be used in this article to indicate the temperature at which the constituent mass of 
quarks vanishes. 
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Pions and the sigma meson are not propogating degrees of freedom in this model, rather they appear 
as poles in the pseudoscalar and scalar channels, respectively, of the effective interaction between quarks 
which is obtained by summing the so-called bubble diagrams. This is the kind of picture that emerges 
for QCD in the limit of a large number of colors N c , so this is the limit in which the NJL model is 
expected to reproduce QCD [6|. This also provides the justification for summing an infinite number of 
diagrams in the gap equation, which contribute at the same order 0( 1), to obtain the nonperturbative 
constituent mass of quarks. 

Previous calculations of the shear viscosity of the NJL model include those that were done dia- 
grammatically mm both at vanishing and finite baryon chemical potentials. In other calculations, 
even though the transport equation was used, the assumption of relaxation time approximation was 
made ® uni in, and the shear viscosity coefficient was calculated in the so-called random phase 
approximation to leading order in N c at vanishing chemical potential |12| . There were other calcula¬ 
tions at vanishing baryo-chemical potential mm which considered the enlarged symmetry group of 
SU(3) V x SU(3 ) a x U(1) B , and that also assumed the relaxation time approximation to be valid. 

There also exist results for the shear viscosity of the quark gluon plasma at leading order of the 
strong coupling constant [2], albeit at vanishing baryon number chemical potential, and at leading 
logarithmic order at finite baryon number chemical potential [1]. These will be compared in Section [I] 
with the results obtained in this paper as well. 

Jeon had shown m in the context of A</> 4 theory that the leading order diagrammatic calculation 
of the viscosity coefficients is equivalent to that done with the Boltzmann equation, and later it was 
shown [16 : , at least in the case of photo production, that this holds true for finite baryon number 
chemical potential, too. It will be assumed in this article that the equivalence between the perturbative 
calculations of the viscosity coefficients diagrammatically and via kinetic theory holds even when the 
baryon number chemical potential is nonzero. 

The aim of this article is twofold: To calculate the shear viscosity of the NJL model at vanishing 
baryon number chemical potential by solving the Boltzmann equation, that is without making the 
relaxation time approximation, and to systematically generalize it to finite baryon number chemical 
potentials. In a leading order in N c calculation, summation of bubble diagrams leads to pion and sigma 
meson exchange in the scalar and pseudoscalar channels, respectively. Above the chiral phase transition 
these appear as resonances, however as the temperature gets closer to the critical temperature the pion 
resonances turn into bound states that may then go on-shell. Although this is innocuous in s-channel 
diagrams, it leads to divergences in t- and u-channel ones. This problem was not present in previous 
calculations utilizing the relaxation time approximation, because an integrated cross-section was used 
there. In order to cut off these divergences one needs to include the width of the pion and sigma 
mesons below and slightly above the critical temperature. However, this width does not appear at 
leading order in N c . We choose to avoid these divergences instead by approximating the interaction 
between quarks with the cross section which is correct to leading order in the coupling constant. 

As ensuing analysis of the NJL model m showed, the interaction between quarks and the mesons 
gets weaker as the temperature approaches the critical temperature from below. Therefore, the neglect 
of bubble diagrams can, at best, be a reasonable approximation as long as the coupling of quarks to 
the effective mesonic degrees of freedom g nqq is less than one. This condition will be satisfied if the 
constituent mass of quarks M is less than the zero-temperature pion decay constant F n ss 95 MeV, 
that is g nqq = M/F\ < 1. This condition determines the lowest temperature at which we will calculate 
the shear viscosity coefficient in this article. 

Dependence of shear viscosity of the NJL model on the temperature will be different depending on 
whether the system is in the chirally symmetric phase or not. Taking into consideration the fact that 
the shear viscosity will be inversely proportional to the square of the coupling constant G[MeV~ 2 ] of 
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Figure 1: Shear viscosity to entropy density as a function of temperature in MeV, shown on the left 
panel are those for the quark number chemical potential to temperature ratios of 0.0, 0.3 and 0.5, 
whereas on the right panel plotted are those for the ratios of 0.6, 0.8 and 1.0. The horizontal dashed 
line on each graph indicates the conjectured lower bound of 77 /s = ( 47 r) _1 , see references [18] . 

the NJL model, in the phase of broken chiral symmetry shear viscosity scales, on dimensional grounds, 
as, 


77 = {G 2 T)~'xf{n/T,M/T), (1) 

whereas for the chirally symmetric phase scaling of the shear viscosity is as follows: 

77 = (G^)- 1 x g^/T) . (2) 

Only two of the three parameters among temperature T, quark chemical potential /t, and M in 
Eq. dTJ are independent. The precise way in which the third parameter depends on the others through 
the gap equation Eq. 0 . and the issue of cutoff dependence of the results overlooked in this section 
to simplify matters will be taken up in Sections [2] and [0 respectively. Section [2] will also give a 
brief description of the NJL model, details of the calculation and how the Boltzmann equation is 
implemented will be explained in Section [3] Section [4] will be devoted to a discussion of our results. 

Our major result in this paper is the calculation of the ratio of the shear viscosity to the entropy 
density for ratios of chemical potential to temperature ranging from /r/T = 0.0 to [i/T = 1.0 as 
shown in Fig. [I] Although it is hard to tell from the figure, the parametric behavior of the ratio of 
shear viscosity to entropy density s changes at the critical values of T c = 190 and T c = 187 MeV for 
/r/T = 0.0 and g/T = 0.3, respectively. For the higher values of the chemical potential to temperature 
ratio, the curves shown in the figure correspond to the chirally symmetric phase only. 


2 Two-Flavor NJL Model 


NJL model is given by the following Lagrangian density: 


C = T 


- 777 .) T + G (’FT) 2 + (Tr^aT)" 


T777 U T 


( 3 ) 


The Lagrangian given above is for Nf = 2, where Nf is the number of quark flavors. Components 
of the matrix vector denoted by a in the second interaction term above are the Pauli matrices. Quarks 
are chosen to be in the fundamental representation of SU(3), then the anti-quarks have to be in the 

2 The values quoted in Fig. [l] are those of the ratio of quark number chemical potential to temperature. 


3 










conjugate representation. Color indices of the quark field will be suppressed throughout the paper as 
neither interaction term in the Lagrangian density depends on them. 

In the chiral limit, the mass matrix is zero m = 0. The last term has to be introduced when the 
chemical potential corresponding to baryon number is nonzero. In order to preserve the symmetries 
of the model it will be chosen to be proportional to the identity matrix (in flavor space), with the 
proportionality constant p being the common chemical potential of both flavors. The relation between 
the chemical potential of the quark species and that for baryon number chemical potential p F is 
M = Atfi/3- 

NJL model in the chiral limit has two phases in the /,i — T plane separated by a second order phase 
transitior|f|. The existence of a nontrivial solution to the gap equation given below determines whether 
the system is in the chirally symmetric high-temperature phase or the low-temperature phase of broken 
chiral symmetrjQ. 


M 


r A d 3 d 

1 - 4GN c N f / - P [1 - n F (T, p) - n F (T, p)} 

Jo A p 


= 0 


(4) 


Here E p = [p 2 + M 2 ) 1 / 2 is the quasiparticle energy and p = |p| denotes the magnitude of the particle’s 
momentum, n F (T,p) = \pJ E p~G/ T ij j s the distribution function of quarks, and h F (T,/x ) = 
n F (T , —p) is that of anti-quarks in thermodynamic equilibrium. Eq. (U) always has the trivial solution 
M = 0 and for a given chemical potential this is the only solution above the critical temperature. 
However, at fixed chemical potential, as the system cools down the ground state that minimizes the 
free energy will lead to a nonzero value of the constituent quark mass, i.e. spontaneous chiral symmetry 
breaking. 

As described in greater detail elsewhere m, an analysis of the NJL model in the chiral limit at 
vanishing temperature and chemical potential determines the value of the coupling constant G and 
the 3-momentum cutoff A in Eq. (J4]) as GA 2 = 2.14 and A = 653 MeV, respectively. Therefore, this 
effective theory can be used for the calculation of shear viscosity for energies below this cutoff scale. 


3 Details of the Calculation 

The cross section for a given particle to scatter off other particles in the plasma scales at@ a ~ G 2 T 2 ~ 
(GT 2 ) 2 T~ 2 . This cross section then leads to a mean free path of l ~ (T 3 (GT 2 ) 2 T -2 ) -1 ~ ( GT 2 )~ 2 T~ 1 . 
Since this mean free path is much larger than the thermal wavelength of quasi-particles in the plasma 
l ~ ( GT 2 )~ 2 T~ 1 T -1 , one is justified in using kinetic theorj@, and more specifically the Boltzmann 
equation to calculate transport coefficients. The above estimate for the mean free path then leads to 
the following estimate for shear viscosity: 77 ~ ^ ~ (GT 2 ) _ 2 T _ 1 T 4 5 ~ (G 2 T) _1 , where £ is the energy 

density of the quark plasma. This estimate was the basis for Eqs. Q and ©. 

We denote the phase space density of both quarks and anti-quarks in the thermal bath by /(x, p). 
In the hydrodynamic limit of small deviations from local thermodyanmic equilibrium one approximates 
this phase space density as a local equilibrium piece / eg (x, p) and a small deviation from that <5/(x, p). 

3 The nature of the transition is very sensitive to the value of the current quark mass, and changes from a second 
order transition to either a first order transition or a cross-over for finite values of this parameter depending on the 
temperature to chemical potential ratio. 

4 The coefficient in the gap equation should actually be N c Nf + however the second term is subleading for large 
N c and is generally dropped. 

5 The estimate given above assumes that T< A, otherwise the scattering cross section will also depend on the cutoff. 

6 In a leading order in N c calculation with N c 1, but GA 2 N C ~ GT 2 N C ~ G( 1) these estimates would have been 
cr ~ (GT 2 N C ) 2 T~ 2 N ^ 1 - l - (GT 2 AT c ) -2 T -1 A/c ~ 0 (N C ) > T _1 , and 77 - ( GT 2 N C )~ 2 T 3 - (9(1). 
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The local equilibrium distribution, of course, has to be taken equal to Uf{T, /i) and np(T, fi) for quarks 
and anti-quarks, respectively. For a divergenceless flow relevant for shear viscosity this phase space 
density can be taken to be independent of time. 


/(x,p) = / e9 (x,p) +<S/(x,p) (5) 

The Boltzmann equation gives the change in the phase space density as a result of collisions, which 
are encapsulated by the collision term C[f] as follows: 

g+v p -V x / = — C[f]. (6) 

The velocity of a particle with momentum p in the above equation is v p = p /E p as usual for a 
relativistic particle. As the LHS of Eq. © is already first order in the gradients of the flow velocity, 
the local equilibrium distribution function should be used there for a leading order calculation. We 
use the invariance of the Boltzmann equation under boosts (more generally Lorentz transformations) 
to evaluate the derivatives on the LHS in a reference frame where the flow velocity is nonzero, and 
specialize to the local rest frame only after this has been done. This gives 


Vp • V x / e ,(x, p) = Pfeq(l ~ feg)^ 

H/p 


-^= + VyUi - -6ijV ■ u 


with the second rank tensor (under rotations) Iij in the above equation defined as 


(7) 


hj{ p)= 



( 8 ) 


Normalization is chosen such that Aj(p)/,;j(k) = P 2 (cos0 p k), where P2(cos0) is the second Legendre 
polynomial. We follow the conventions of [2] in the definition and normalization of the above tensor. 
The departure from equilibrium at linearized order can naturally be parametrized as follows: 


<5T(p)=/3 2 /e>) 


[V6 


ViUi 


Vy Uj ’ U 


Xl 


(9) 


This leads to the simpler equation given below. The expression for the linearized collision operator C 
appearing in the equation below will be given later. 


<S?j(p) = {CXijTi P) (10) 

where the source term 5“ is defined to be 

L J 

S ? 3 = - T ^/e>) [1 - /e»] Iij (P) (11) 

and due to the rotational invariance of the collision operator, once the tensorial structure of \ij is 
peeled off it can only be a function of the magnitude of momentum in the local rest frame. 


Xij = Iij(p)x a (p) 


( 12 ) 


It’s convenient to define the following inner product 

(= j [1 - feq(p)\ 


(13) 
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where we used Iij(p)Iij{p) = 1. Expanding the solution in terms of basis functions given below 

K 

X a {p)=^2x a rn<t><,rn){v) (14) 

m =1 

gives the vector representation of the source and the matrix representation of the linearized collision 
operator in this basis as follows. 


(Xij i &ij ) — 

{Xij ) Cxij ) = 

In Eq. (fT4j). index a denotes the type of particle, in our case either quark or anti-quark, for which this 
is the deviation from equilibrium and the lower index stands for the coefficient of the basis function m 
in the expansion of the deviation from equilibrium. Here is the explicit expression that we postponed 
to give for the linearized collision operator C^ n appearing above. 

C a ± = lY, f f^E^f^E^ll-f^E^ll-f^E^lM ^ 2 

8 a,6,c,d^P’ k ’P'> k ' 

X [r { mMV>) + ~ ^mM^) 

■ + 0(„)^(k) - 4*) W) - 

x (2 tt) 4 <5 (4) (P + K-P' - K') (16) 


V Y° S a 

/ j A.m^m 
a,m 

(15) 

a,m b,n 


The above integrals can be simplified considerably, details of the simplified collision term are given in 
Appendix |Aj One basis that has the right asymptotic behavior and converges with a sufficiently small 
number of basis elements is provided by the following functions. 


^(m) ( P ) 


v{vlT) m - x 

(l+p/Tf- 2 


(17) 


with m = 1,2In our case, when the cutoff is neglected choosing five basis elements, that is 
K = 5, leads to one part in a thousand accuracy, whereas inclusion of the cutoff severely reduces this 
to slightly above 1.5%. Shear viscosity is then given in terms of the source vector and the inverse of 
the collision matrix in this basis as follows: 

= = (18) 

A few comments are in order about the matrix elements used in Eq. (HU), in the collision term of 
the Boltzmann equation. Since we have worked to leading order in N c and treated GA 2 N C ~ 0(1) 
in the gap equation, consistency would require that we include in the matrix elements all scattering 
diagrams to leading order in N c , that is the so-called bubble diagrams, for a given process as well. 
These diagrams, when summed, would then lead to the effective interaction of a pion being exchanged 
in the pseudo-scalar channel, and a sigma meson being exchanged in the scalar channel. Although 
this is the correct procedure to follow, we will ignore the pion (or the sigma meson) pole contribution, 
and approximate the quark and anti-quark interactions by the tree level term read off from the NJL 
model lagrangian. We are going to compare our results for the zero baryon chemical potential case to 
previous work in Section [5] to see how this approximation affects our results. 

It is in principle possible to include the pion pole contribution to the scattering matrix elements 
in the calculation of the shear viscosity of the NJL model, however it should be noted that this is 


6 







M a c b d 

\M“ b \ 2 /G 2 

qq qq 

M qq (P, K;P', K') 

256V 2 (t 2 + u 2 — tu + 6M 4 ) 

99 ^ qq 

M qq -(P , K; P', K') = M qq (P', K'- P, K) 

256 N 2 (t 2 + u 2 —tu + 6M 4 ) 

99 ^ 99 

M qq (P , K; P', K') = M qq (P, —K'; P', -K) 

256V 2 (t 2 + s 2 — ts + 6M 4 ) 


Table 1: Charge conjugation invariance forces the matrix elements for scattering of quarks off quarks 
and of antiquarks off antiquarks to be the same, whereas crossing symmetry relates the matrix elements 
for scattering of quarks off quarks to scattering of quarks off antiquarks. The last column gives 
the absolute square of the matrix elements divided by the square of the coupling constant for these 
processes. The square of the matrix elements are summed, but not averaged, over the spins, colors and 
flavors of the participants. In particular, the matrix elements squred scale wih V 2 because an average 
over the colors has not been taken. Of course, the number of colors will be set equal to 3. 

nontrivial due to the divergences it causes in the collision term of the Boltzmann equation. We work 
in the chiral limit throughout the paper, and in this limit the pion is exactly massless. Therefore, the 
divergences in the s-channel pion exchange is harmless due to shrinking of the phase space, and a bare 
pion propogator is sufficient for s-channel exchange. This parallels s-channel gluon exchange in a quark 
gluon plasma calculation of shear viscosity where it suffices to use a bare gluon propogator. However, 
t- and u-channel pion exchange lead to divergences unless the imaginary part of the pion propogator, 
that is generated by scatterings of the pion off quarks and anti-quarks in the plasma, is included in 
those matrix elements. This is expected to arise at order 0{N ~ 2 ) for the cross-sections, that is at 
next-to-leading order in N c . 

As for the masses of the species in the matrix elements and in the equilibrium Fermi-Dirac distri¬ 
bution functions, these are taken to be zero in the chirally symmetric phase, and non-zero and equal 
to the value obtained from Eq. Q for a given temperature and chemical potential in the calculation of 
shear viscosity explained in this paper. For our parameter choice, the critical temperatures for chiral 
symmetry breaking at /i/T = 0.0 and /r/T = 0.3 are T c = 190 and T c = 187 MeV, respectively. The 
critical temperature is below 185 MeV for all the higher ratios of the chemical potential to temperature 
ratio, so the quasi-particles are taken to be massless in those cases. 

It should be emphasized that the masses of the quasiparticles in the NJL model are nonperturbative 
masses, and thus they are included in the low temperature phase where they are nonzero. The treatment 
of masses in this model should be contrasted with how the thermal masses of species were handled in 
a leading order in the coupling constant calculation of the shear viscosity of the quark gluon plasma 
m 2 h ~ 9sT 2 and of scalar theory with quartic self-couplings m 2 h ~ AT 2 , where g s and A are the 
coupling constants of QCD and A(/> 4 theory, respectively. In those calculations, all on-shell particles 
were taken to be exactly massless, and the perturbative effects induced by the thermal masses rightly 
belonged among the corrections higher order in the coupling constant the leading order result receives. 

Charge conjugation invariance of the kinetic and interaction terms in the NJL model, that is barring 
the chemical potential term, implies that the matrix elements for the scattering of quarks off quarks 
and of anti-quarks off anti-quarks are the same. Moreover, the matrix elements for the scattering of 
quarks off anti-quarks is related to the above by crossing symmetry. Explicit expressions, in terms of 
the Mandelstam variables, for the squares of these matrix elements as they appear in the Boltzmann 
equation, that is summed but not averaged over the spins, flavors, and colors of the participants, are 
given in the Table 1. 
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Figure 2: The dimensionless product of shear viscosity, coupling constant squared and the temperature 
as a function of temperature in MeV. The horizontal lines in the plot show the constant values of this 
product when the cutoff is taken to infinity for the chemical potential to temperature ratios of 0 . 6 , 0.8 
and 1.0. The scattered points, on the other hand, are data for the same ratios of the chemical potential 
to temperature when a momentum cutoff is used in the Boltzmann equation. 


4 Results &; Discussion 


In order to understand the manner ratio of shear viscosity to temperature cubed depends on tem¬ 
perature, we begin by analyzing the dimensionless product r?G 2 T, this product is expected to only 
depend on the ratio of chemical potential to temperature in the chirally symmetric phase, and should 
therefore be constant for each value of this ratio. This is clearly the case when the 3-momentum cutoff 
is taken to infinity in the Boltzmann equation as can be seen in Figure [2] for various values of the 
chemical potential to temperature ratio. Thus, the ratio of shear viscosity to temperature cubed scales 
as r]/T 3 ~ T -4 in the high temperature phase. The reason the ratio of shear viscosity to the cube of 
the temperature varies with temperature at all in the massless phase is because the effective coupling 
constant GT 2 = GA 2 (T/A ) 2 = 2.14(T/A) 2 , which is what is relevant for scattering cross-sections, 
increases with increasing temperature. 

We evaluated the shear viscosity coefficient in the phase of broken chiral symmetry only for the 
chemical potential to temperature ratios of 0.0 and 0.3, as the transition temperature in all the other 
cases falls below 185 MeV. This is the temperature for the vanishing chemical potential case at which 
the pion-quark coupling constant g^qq given by the ratio of the constituent mass of quarks to the zero 
temperature pion decay constant F n is of order one g^ qq = M/F v ~ 0( 1), and the ’’fine structure 
constant 1 ' is negligible a = g 2 qq /(ATr). Our calculation can, at best, be expected to be a good approxi¬ 
mation down to these temperatures. Since this temperature is very close to the transition temperatures 
of 190 and 187 MeV for the vanishing chemical potential and for the chemical potential to temperature 
ratio of 0.3, respectively, the dimensionless product of the shear viscosity, coupling constant squared 
and the temperature r)G 2 T can be Taylor expanded around its value in the chirally symmetric phase 
tiG 2 Tt<t c ~ vG 2 T c — a ■ ( T c — T) with a positive constant a > 0. We found that the values of this 
dimensionless product were fitted very well with a line in the low temperature phase, the slope of this 
best fit is the above constant a. Ratio of shear viscosity to the cube of the temperature, unlike the 
dimensionless product, increases as the temperature is lowered in the massive phase. However, the 
scaling with temperature changes to T f^ Tc + 7 ^ 3 . The change in the scaling of the ratio of 

shear viscosity to temperature cubed in the two phases is apparent as there is a discontinuity in the 
slope of the lowest two curves in the left panel of Fig. [3] at the transition points of 187 and 190 MeV. 
These discontinuities persist in the lowest two curves in the left panel of Fig. [TJ too, though they are 







not as prominent. 

Shear viscosity of the NJL model at vanishing chemical potential was calculated previously by 
Zhuang et al. (T2J. What we found for the ratio of shear viscosity to temperature cubed in the 
low temperature phase of broken symmetry agrees with their result exactly. Although they use the 
same matrix elements in the massive phase as we do, this could not be expected a priori given the 
relaxation time approximation they made in their calculation. On the other hand, the dependence 
of this ratio on the temperature is quite different for the two calculations in the chirally symmetric 
high temperature phase. Shear viscosity scaled with the third power of the temperature decreases 
with increasing temperature in both calculations, yet the scalings are different. In this paper, we find 
this ratio to decrease with the fourth power of the temperature r)/T 3 ~ T~ 4 , whereas the figure in the 
paper of Zhuang et al. indicates a less steep decrease with increasing temperature than ours. Thus, our 
calculation underestimates this ratio by a factor of nearly 3.5 compared to theirs at the temperature 
of 300 MeV. 

Notwithstanding the fact that the ratio of shear viscosity to entropy density was calculated previ¬ 
ously, it should be emphasized that this ratio calculated with the NJL model accounts for the almost 
saturation of the lower bound of this value observed in experiments conducted at RHIC much better 
than perturbative QCD calculations do. For example, the leading order in the strong coupling constant 
calculation of Arnold et al. [2: found the viscosity coefficient to be 77 = /(tod/F) where mn is the 
Debye mass. The value of the ratio of Debye mass to temperature for two flavors considered in this 
paper is approximately mo/T ps 2.24 for the phenomenological value of a s cs 0.3. Figure 2 in the 
aforementioned reference gives for the value of this function fimo/T ps 2.24) a; 150 so that the ratio of 
shear viscosity to entropy density comes out to be r\/s ps 2.60 which is nearly 3 times the highest value 
of this ratio obtained in this paper for the vanishing chemical potential case. On the other hand, the 
lowest value of the shear viscosity to entropy density is attained at the upper end of the temperature 
range, and is 77 /s ps 0.13 which is only 60% greater than the conjectured lower bound |T8]. The scalings 
of the ratio of shear viscosity to the cube of the temperature is also different in the two calculations, 
high temperature perturbation theory predicts this ratio to be independent of temperature, save for 
the logarithmic running of the coupling constant with temperature, whereas in this article we found 
this ratio to decrease with the fourth power of the temperature r]/T 3 ~ T -4 . Finally, the ratio of shear 
viscosity to entropy density increases monotonically as the ratio of chemical potential to temperature 
is increased from n/T = 0.0 to fi/T = 1.0. 

If the temperatures at which we evaluated the shear viscosity coefficient had been much smaller 
than the 3-momentum cutoff TC A that appears in the gap equation Eq.(|4]), and which defines the 
theory, the precise value of the cutoff would have been immaterial. This is because the equilibrium 
distribution functions in the collision term of the Boltzmann equation decay exponentially for energies 
much larger than the temperature E p T and particles with such high momenta do not contribute 
to the integrals. However, for the temperatures of 185 to 300 MeV at which we are interested in 
calculating the shear viscosity coefficient one has to assess how sensitive one’s results are to the value 
of this cutoff. This is exactly what the scattered points in Figure [2] indicate. As can be seen from 
the figure, at lower temperatures imposing a 3-momentum cutoff tends to increase the dimensionless 
product, and thus the viscosity coefficient itself, whereas at higher temperatures it tends to reduce the 
value of the shear viscosity to almost half the value it has without the 3-momentum cutoff. 

In this article we calculated the leading order in the coupling constant shear viscosity coefficient of 
a plasma of quarks and anti-quarks interacting according to the NJL model at finite temperature above 
and slightly below the chiral phase transition temperature and at nonzero baryon number chemical 
potential. In the calculation the ratio of chemical potential to temperature was not greater than one, i.e. 

7 m 2 D = yr( 3 + N f /2)a s T 2 
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Figure 3: Ratio of shear viscosity to temperature cubed as a function of temperature in MeV. The 
plot on the left panel shows the curves for the quark chemical potential to temperature ratios of 0 . 0 , 
0.3 and 0.5, whereas on the right panel displayed are those for 0.6, 0.8 and 1.0. Note the discontinuity 
near 190 MeV in the slope of the lower two curves shown on the left panel. 

n/T < 1. We did not sum all scattering diagrams to leading order in the number of colors 0(N~ 1 ), 
and instead evaluated the matrix elements for scattering by the tree level scalar and pseudoscalar 
interaction terms in the Lagrangian, that is to leading order in the coupling constant 0(G 2 ). In 
order to obtain the shear viscosity coefficient we solved the linearized Boltzmann equation. We found 
that the ratio of shear viscosity to entropy density decreased with increasing temperature, scaling like 
77 /s ~ T ~ 4 in the high temperature phase. The value of this ratio was found to be roughly 1.5 times 
the conjectured lower bound of (47 t) _ 1 at the temperature of 300 MeV for the vanishing chemical 
potential case, and was not in any of the cases considered greater than 13 times this bound. Therefore, 
shear viscosity over entropy density calculated in this model reproduced the low value observed in 
experiments far better than perturbative calculations did. Shear viscosity did not change much as 
the ratio of chemical potential to temperature was varied, however it was a monotonically increasing 
function of this variable. 

The vanishing width of the pion as the critical temperature is approached from above prevented 
us from using all cross sections for scattering that are correct to order 0(N~ 1 ) in the Boltzmann 
equation. As even stable particles acquire a width in a thermal medium, it might be possible to 
correctly incorporate this width of the pion in the calculation of the shear viscosity by including 
diagrams of 0(N~ 2 ). This nontrivial extension of the current calculation of shear viscosity is left for 
future work. 


A Integration Variables and Limits 

In order to simplify Eq. (TTflt we follow the steps outlined in the appendix of Arnold, et al. [2], but 
generalize their treatment to take into account the finite mass of quarks. For terms that are proportional 
to t = {P — P') 2 the spatial delta function can be used to perform the k' integration and the p' 
integration can then be shifted into one over q = p' p. This reduces Eq. m to the following form 

o3 r °° /'I />27T I a /iab\2 

iXij,C 2 ^ 2 Xij) = 77 TJ V / q 2 dqp 2 dpk 2 dk d cos6 pq d cos6 kq # g* x 

(4tt J abcd do J -1 JO E p E k E p 'E k ' 

S(E P + E k - E p> - E k .) r eq {E p ) f b eq (E k ) [1 - f c eq (E p ,)] [1 - f d eq {E k ,)\ x 

[Xij (P) + Xij (k) - Xij (pO - xtj (k')] 2 (19) 
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At this point a new integration variable w, which corresponds physically to the energy transferred 
in the collision, is introduced following Baym et al. 20j to do the angular integrals as follows. 

/ +oo 

dui S (uj + E p — Epi) 6 ( w — E k + E k >) (20) 

-OO 

Using the properties of the delta function one finds that the above delta functions can be written 
instead as 


6 {uj + Ep — Ep /) 


5 (w - E k + E k ') 


Ep' 

pq 

kq 


cos 9p Q - 

q 

cos 6 kq - — 



2S) e( "’ + E " ) 

i) 0(Et - w> 


( 21 ) 


with 0 denoting the unit step function. Carrying out these delta function integrals reduces Eq. m 
further to the form below 


(Xtf.C^Xy) 


o 3 _ /*oo r w + r°° r 00 /*27t 

71^6 J2 dq dw dp dk ^ \ M cd\ 2 
(47T) a bcci ^ Jw- Jp~ J k- Jo 

feq(Ep) f b eq (E k ) [1 - f^Ep')} [1 - f« q (E k ,)\ X 

[Xij( p) + Xij{ k) - Xij( p') - xfj( k ')] 2 



( 22 ) 


where the upper and lower limits of the w integral are w± = ± ^ \Jq 1 + M 2 — M ^, the lower limit of the 

p integral is p- = (q — w\Jl — 4 M 2 /t^j /2, and that of the k integral is = (q + w\J 1 — 4 M 2 /t^j /2. 
In the above equation, E p ' = E p + w and E k ' = E k — w. The Mandelstam variable t is as usual 
t = w 2 — q 2 . Imposing a 3-momentum cutoff in the Boltzmann equation is implemented by setting the 
upper limit on the q integral in Eq. (1221) above to twice the value of the 3-momentum cutoff A. 

In the t-channel parametrization, the other Mandelstam variables s and u are related to t as given 
below. 


s — 2 M 2 — ^2 [(E P + Ep')(E k + E^) + q 2 ] - 

x^/[{-t)(AE p Ep' + t) - Aq 2 M 2 \ [(-t)(4E k E k . + t) - 4 q 2 M 2 ] , (23) 

u = AM 2 -s-t , (24) 


where M is the constituent mass of the u and d quarks. Finally, one needs the cosine of the angles 
between various momenta and these are given below. 


cos 0 P q = 



cos 0pp' — 
cos Opk' — 


EpEp' - m 2 
pp' 

E p E k f - m 2 
pk' 


t 

2 pq ’ 


t 

2 pp' 


u 

’ 



w 

q 


t 

2 kq 


cos 0 kk ' — 


E k E k ' - m 2 
kk' 


t 

2kE ’ 


cos0p' k — 


E p rE k - m 2 
p'k 


u 

2 p'k 


(25) 

(26) 
(27) 


The neccessary equations for the chirally symmetric phase in which the constituent mass of quarks 
vanishes can be obtained from those above by simply setting M = 0. There are equations similar to 
the ones above for the s-channel parametrization, however those will not be reproduced here. 


11 

















References 

[1] P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0011, 001 (2000). 

[2] P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0305, 051 (2003). 

[3] P. B. Arnold, Q. Dogan and G. D. Moore, Phys. Rev. D 74, 085021 (2006). 

[4] J. W. Chen, Y. F. Liu, Y. K. Song and Q. Wang, Phys. Rev. D 87, 036002 (2013) [arXiv:1212.5308 
[hep-ph]]. 

[5] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); Phys. Rev. 124, 246 (1961). 

[6] A. H. Blin, B. Hiller and M. Scliaden, Z. Phys. A 331, 75 (1988); E. Witten, Nucl. Phys. B 160, 
57 (1979). 

[7] Y. Hidaka and T. Kunihiro, J. Phys. Conf. Ser. 270, 012050 (2011). 

[8] R. Lang and W. Weise, Eur. Phys. J. A 50, 63 (2014). 

[9] S. Ghosh, A. Lahiri, S. Majumder, R. Ray and S. K. Ghosh, Phys. Rev. C 88, 068201 (2013). 

[10] C. Sasaki and K. Redlich, Nucl. Phys. A 832, 62 (2010). 

[11] M. Buballa, K. Heckmann and J. Wambach, Prog. Part. Nucl. Phys. 67, 348 (2012). 

[12] P. Zhuang, J. Hufner, S. P. Klevansky and L. Neise, Phys. Rev. D 51, 3728 (1995). 

[13] R. Marty, E. Bratkovskaya, W. Cassing, J. Aichelin and H. Berrehrah, Phys. Rev. C 88, 045204 
(2013). 

[14] P. Rehberg, S. P. Klevansky and J. Hufner, Nucl. Phys. A 608, 356 (1996) hep-ph/9607263|. 

[15] S. Jeon, Phys. Rev. D 52, 3591 (1995); S. Jeon and L. G. Yaffe, Phys. Rev. D 53, 5799 (1996). 

[16] H. Gervais and S. Jeon, Phys. Rev. C 86, 034904 (2012). 

[17] V. Bernard, U. G. Meissner and I. Zahed, Phys. Rev. D 36, 819 (1987); T. Hatsuda and T. Ku¬ 
nihiro, Phys. Lett. B 185, 304 (1987). 

[18] P. Kovtun, D. T. Son and A. O. Starinets, JHEP 0310, 064 (2003); Phys. Rev. Lett. 94, 111601 
(2005). 

[19] S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

[20] G. Baym, H. Monien, C. J. Pethick and D. G. Ravenhall, Phys. Rev. Lett. 64, 1867 (1990). 


12 


